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ABSTRACT 


Modern telecommunication scenarios, especially those related to 5G, pose significant 
challenges to electromagnetic simulation tools. The scenarios are extremely large in terms of 
wavelengths and are geometrically complicated. Even ray-tracing methods can easily require 
too much simulation time, as too many rays need to be traced. This presentation was given as a 
keynote speech at MOMAG 2018. It compares several electromagnetic simulation methods for 
indoor, urban and rural propagation, and discusses the trade-offs in accuracy and simulation 
time. 


Index terms: Electromagnetic simulation, propagation, 5G, ray tracing 
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EVOLUTION OF CELLULAR GENERATIONS 
1G: mostly analog 2G: digital; Voice + texting; FDMA and TDMA 





4G: Higher data rates 
with MIMO and OFDM 





CHALLENGES FOR 5G 


5G Service and Scenario Requirements 

* Massive growth in mobile data demand (1000x capacity) 

* Higher data rates per user (10x) 

* Massive growth of connected devices (10-100x) 

* Lower latency for real-time connections (5x) 

* Higher energy efficiency for longer battery life (100x) 

* New use cases like connected cars, machine-to-machine (loT) 


Key Solutions 


e Increased spectrum with trend to higher frequencies 
e 700 MHz for basic coverage 
e 3.5 GHz for high data rate services & capacity 
¢ 26 GHz for fiber-like data rates & capacity hotspots 


ο Ultra-dense networks 


e Massive MIMO antennas for beamforming & spatial 
multiplexing 


Multi-User UHD 
Telepresence 
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Source: Huawei 5G Technology Vision 


5G: MAIN CHARACTERISTICS 
* Higher frequencies, e.g. 26 GHz and above 
* Massive MIMO 


e Beams that follow the user 
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CEM SOLVER TECHNOLOGIES 
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COMPLEXITY OF MATERIALS 
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THE DEVICE 
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DEVICE: FROM ANTENNA INTEGRATION TO ANTENNA PLACEMENT 











Ε-[Ι8ία 0.9 A simulation time: 43min Е-Неа 26 4ΗΠΖ ΠΠ ο time: 56min 


Previously, extremely detailed CAD geometry was cumbersome 
* Over-discretize the FDTD mesh to resolve geometric detail 


Now default meshing is « 1mm, most detail is inherently captured 
* Despite the electrical size at 26 GHz the integrated antenna simulation can be run in « 1hr 


PCB: part of antenna at low frequency — large ground plane at 26 GHz 
This phone with small array will be used later in 5G network planning example. 
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THE BASE STATION 








Gain [dBi] 


ARRAY DESIGN - GAIN & BEAM STEERING 


4x4 array array 8x8 array 16x16 












Gain [dBi] 


Theta [deg] 


gain for 4x4, 8x8, 16x16 array configurations 


This example array will be used later 
in 5G network planning example. 
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8x8 array 
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Theta [deg] 
beam steering for the 8x8 array 


RADIO CHANNEL AND COVERAGE ANALYSIS 
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Hong Kong, 
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EXAMPLE: COMBINED URBAN+INDOOR 


Detailed multi-floor 
indoor building 
in urban surrounding 
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EXAMPLE LARGE COMPUTATIONAL DOMAIN AND COMPLEXITY 






Multi-floor 
building 
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FEKO - WINPROP INTERACTION 


Not feasible with SBR 
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3D pattern for three-sector antenna computed in FEKO | τὰ, | 4 
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eet 


We will discuss key WinProp methods: 
- Standard Ray Tracing; Intelligent Ray Tracing 
- Dominant-Path Model 


Εξ „йы J LS = 


Urban radio coverage considering this antenna computed in WinProp 





PROPAGATION SIMULATION METHODS 
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¢ Possibility to define areas with non-deterministic objects (polygonal cylinders) 
* Propagation paths inside these areas get higher attenuation. 





PN 


Without furniture Deterministic Non-deterministic 
and persons modeling modeling 
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INDOOR PROPAGATION MODELING 


Propagation Models 


* One-Slope Model 
e Only distance dependency 
e Path loss exponent is dominant parameter 


e COST 231 Multi-Wall-Model 
e Only direct ray between Tx and Rx 
* Transmission loss for each wall COST 231 Multi Wall Ray Tracing 





e Ray [racing 
e 3D Ray Tracing (IRT, with preprocessing) 
e 3D Ray Tracing (SRT, without preprocessing) 


e Dominant Path Model 
. 3D (Multiple floors) 


Dominant Path 
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MULTI-WALL MODEL 


Principle of model: 

e Based on direct ray between transmitter and receiver 
e Free space loss (one slope) with additional loss due to transmission of walls 
e Individual material properties of each wall considered 


Field Strength 
(dB uV/m) 
120.00 
Receiver 110.00 
100.00 
90.00 
80.00 


70.00 


Computation of field strength: 


Transmitter 60.00 


50.00 


40.00 


30.00 





20.00 
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RAY TRACING 


Multipath propagation considered 


Dominant effects: diffraction, reflection and 
transmission (penetration) 


Full 3D approach 


* Uncorrelated or correlated superposition of 
contributions (rays) 


e [wo sub-modes: 


e 3D IRT: with pre-processing of building data 
and accelerated predictions 


e 3D SRT: without pre-processing 
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PRE-PROCESSING FOR INTELLIGENT RAY TRACING 


* Considerations to accelerate the time consuming 
process of path finding: 


* Visibility relations between objects are determined 
once and used for any transmitter position. 


* Deterministic modelling generates a large number 
of rays, but only few of them deliver most of the 
energy. 


e Adjacent receiver pixels are reached by rays with 
only slightly different paths 
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DOMINANT PATH MODEL 





e Dominant Path (single path or group of paths) 


* Unlimited number of interactions (changes of 
orientation) along the path 


* Parameters of path determined (e.g., length, number 
of interactions, angles...) and used to compute path 
loss with semi-deterministic equations 


* Optional consideration of wave guiding possible (wave 
guiding factor, based on reflection loss of walls) 


* No pre-processing required 











* Accurate and fast 
* Auto-calibration available 
e Robust against errors in vector building database 





empirical model ray-optical model dominant path model y 
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INDOOR RESULT COMPARISON 
Institute for Radio Frequency Technology, University of Stuttgart, Germany 


Field Strength 
[dB V/m] 
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Multi-wall Model 3D Ray Tracing Dominant Path Model y 


DEVIATION FROM MEASURED RESULTS 





Multi-Wall Model 3D Ray Tracing Dominant Path Model y 


© 2018 Altair Engineering, Inc. Proprietary and Confidential. All rights reserved. 


URBAN DATABASES: MORE BUILDINGS, LESS DETAIL PER BUILDING 


e Buildings are extruded polygons 
e > Vertical walls and flat roofs 


e Individual material properties of 
building surfaces 


° Topography can be considered 
optionally 





Example: New York 
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EMPIRICAL METHOD: COST 231 WALFISCH-IKEGAMI 


Power 
[dBm] 


* Model accepted by ITU-R 

e Evaluates building profile between 
transmitter and receiver (vertical plane) 

e Includes additional losses due to building 
data 

e Reasonable results for Tx above rooftops 
For Tx below rooftops limited accuracy (no 
wave guiding) 

e No multipath propagation considered 
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URBAN RESULT COMPARISON 


Helsinki, Finland 
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Empirical 3D Ray Tracing Urban Dominant Path 7 Y 
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Urban Dominant Path 7 Y 


3D Ray Tracin 


Clutter database specifies land usage. 


* Example: Germany 
e 12 classes 
ο 50m resolution 





d Confidential. All rights reserved. 


Clutter Morpho 
[Class | 


water 
ll cpenarea 
Bl parks 
E] forest 
М орегіпиграп 
| industrial 
B ua 
лсе 
> residential 
Ш mearurban 
Щ derseurban 
Мол: 
Ш unclassified 


F 


w 


© 2018 Altair Engineering, Inc. Proprietary and Confidential. All rights reserved. 


TOPOGRAPHY AND CLUTTER DATA 


e Individual class assigned to each pixel} 
Class ID with individual properties Ё 
e Frequency dependent attenuation 
e Clutter heights 

e Clutter clearance 

* Electrical properties of ground 


Clutter/Morpho 
[Class] 


|. Transportation 
Т Airport 

ΚΙ Core Urban 

j Dense Urban 

В urban 

ll Commerc. Indust 
I) Heavy Indust. Re 
~ Light Dens. Res. 
I) Coniferous For. 
I) Deciduous For. 


E Open 


B water 





Example: Detroit, USA 
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WAVE PROPAGATION MODELS FOR RURAL SCENARIOS 


¢ Empirical models (e.g. Hata-Okumura) 


° Only consideration of effective antenna height 
(no topography between Тх and Rx) 2D Vertical plane 


3D Paths 

















* Considering additional losses due to clutter data 


° Semi-Empirical models (e.g. Tvvo-Ray plus 
Rnife-Edge diffraction) 


e Including terrain profile between Tx and Rx 


* Considering additional losses due to diffraction 














e Deterministic models (e.g. Ray Tracing) 
e Considering topography 





e Evaluating additional obstacles ΠΤΙ ο Rx2] 
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TOPOGRAPHY AND CLUTTER DATA 


Propagation Models: Two-Ray & Knife-Edge Diffraction 
* Computation of direct and ground reflected ray 
e Additional diffraction loss in shadowed areas (frequency dependent) 
e Topography between Tx and Rx considered 
e (e.g. shadowing due to hills, etc.) 
° Possible evaluation of Fresnel zone 





ü 10000 20000 30000 40000 50000 


Distance [m] 
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RURAL EVALUATION 


Area around Grab/Murrhardt, Germany 


Topography 
[m] 
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3D view of the database (z-axis scaled with factor 5) 
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RURAL RESULT COMPARISON 





Field Strength 
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Hata-Okumura Model with | Rural 
Knife-Edge-Diffraction Extension Dominant Path Model 
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RURAL DEVIATION FROM MEASURED RESULTS 


Difference 
[dB] 
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Hata-Okumura Model with Rural 
Knife-Edge-Diffraction Extension Dominant Path Model 
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VVAVE PROPAGATION CONSIDERATIONS 
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RAY TRACING 


Multipath Propagation 
e Multiple propagation paths between Tx and Rx 


e Shadowing, reflection, diffraction, scattering 

e Different delays and attenuations 

e Destructive and constructive interference 

e Depending on frequency 

e Various bands of interest for 5G: 700 MHz, 3.5 GHz, 26 GHz, ... 


Superposition of Multiple Paths 
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No line of sight (Rayleigh fading) Line of sight (Rice fading) 


WAVE PROPAGATION » 6 GHZ 


Coverage for Tx Below Rooftop Level (as in 5G) 


Multipath situation 
e Multiple reflections 
e Wave guiding in street canyon 


e Few rays over the rooftops (diffraction) 
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WAVE PROPAGATION > 6 GHZ 


Impact on interactions due to higher frequency 
¢ Transmission 


cy 


Altair 


* Penetration of walls hardly feasible 

e LOS and OLOS regions will dominate (impact of street grid) 
e Reflection 

e Specular paths will dominate (besides direct path) 

e For LOS and OLOS regions 
e Diffraction 





e Highly attenuated for higher frequencies as diffraction coefficient ~ 1/sart(frequency) 
° Will more and more disappear for frequencies > 26 GHz 
e Scattering 
e Roughness becomes large for most surfaces (due to small wavelength) = diffuse scattering 
e 5G transmission will use highly directive antennas on both ends » scattering difficult to be used for reliable connection 


> 5G will need higher density of base stations PN 
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WAVE PROPAGATION > 6 GHZ 


Atmospheric absorption and rain attenuation at mm-wave frequencies 
¢ Additional attenuation tolerable for cell sizes on the order of 200m 
« Atmospheric loss « 0.1 dB/km at 30 GHz, but 20 dB/km at 60 GHz due to oxygen absorption 





e Rain attenuation limited for frequency bands around 26 GHz and 28 GHz 100 ————=—————— 
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oG RADIO CHANNEL 
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9G RADIO CHANNEL 


Massive MIMO antenna arrays 

e Arrays with hundreds of antennas for 
separating tens of users in same radio 
resources (time/frequency) 
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e At mm waves large arrays are compact 


e Combination of 


e Beamforming | | сы 
Massive Arrays of 128-1024 Active Antenna Elements 


UEZ; 16 Ant BF with 8 x 2 MIMO 


* Spatial Multiplexing j R. Ἐν πο. 

Fiber. Fiber. Fiber 
° Relevant channel statistics Ë 
Increase SINR and Canet for Each U = 
* Delay spread Le, UET: 32 Ant BF with 16 x 2 MIMO - = 


e Azimuth angular spread 
both for BS and MS 


e Elevation angular spread Source: Microwave Journal March 2017 





» Evaluation of cumulative 
distribution function (CDF) 
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* Used for prediction of radio channel 
in time, frequency and spatial domain 


* Field strength 
* Propagation delays 
* Angles at Tx and Rx 
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Cummulated Probabüty 
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TYPICAL RESULTS doo ——— x mE p 















Cumulative Channel Statistics 


Cummulated Probability 
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e Computed for individual cells 
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* Consideration of omni BS antenna Dela 
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Cell areas given by best server map 
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* Delay spread at 26 GHz 


Azimuth angular spreads at 26 GHz 
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e Massive MIMO antenna arrays used to | 


* [ransmit different signals to 


different users simultaneously 7 


in same frequency band ке al Z, | ον. us 


° Increase Rx power levels and 
SNIR for dedicated user 


e Reduce interference for others 


e Example shows separation 
of two different users 
* 4x4 antenna elements on BS side 
(upper figures) 
* 16x16 antenna elements on BS side 
(lower figures) 
* [rade-off between beam width and N 
MIMO spatial multiplexing within beam | № 





© 2018 Altair Engineering, Inc. Proprietary and Confidential. All rights reserved. 


PATH LOSS 4-3 — 


-90.00 
-100.00 


5G Radio Channel Analysis | μα И de M o L -110.00 


-120.00 


* Comparison of simulated path loss at 28 GHz P e E, 22| and 2.9 GHz | 


-150.00 


* New York city scenario ala RET E -160.00 


-170.00 


¢ WinProp 3D ray tracing model 
e Base Stations at street intersections 





e 3 different BS locations analyzed 

e Areas marked in black rectangles evaluated in below diagram 

¢ Path loss over BS = MS distance gives much smaller range for 2.9 GHz 
e Wide range of path loss for same/similar distances at 28 GHz 





T. 
U 
e Hence simple path loss model not sufficient 2 
x . © - 
but 3D ray-optical model required Ξ 
σα - 
e Alternative to measurements is the usage of 3D ray tracing engine > 
Source: Qualcomm Z. Zhang et al.: Coverage and Channel Characteristics 
of Millimeter Wave Band Using Ray Tracing, IEEE ICC 2015 Tx-Rx Separation Distance (m) 


© 2018 Altair Engineering, Inc. Proprietary and Confidential. All rights reserved. 


2nd beam 


9G RADIO CHANNEL: BEAM DIVERSITY 


5G Radio Channel Analysis II 
* Beam diversity at 28 GHz 
* New York city scenario 
e WinProp 3D ray tracing model 
* Base Stations at street intersections 


° Three different Base Station locations 
analyzed 


e Statistics on number of paths between 
Base Station and Mobile Station 
> typically (when CDF=0.5) 
there are 2-4 clusters (see figure 1) 


Empirical CDF of Number of Clusters 





10 12 14 


6 5 
Number of Clusters 


Source: Qualcomm Z. Zhang et al.: Coverage and Channel Characteristics 


of Millimeter Wave Band Using Ray Tracing, IEEE ICC 2015 


Figure 1 C 
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BEAM DIVERSITY, CONTINUED 


e |f dominant direction is blocked alternative path required 
* Power fraction of the second strongest path on average 7 dB weaker (Fig. 2) 
¢ Azimuthal separation between two strongest paths on average about 20? (Fig. 3) 


¢ WinProp simulation results in agreement with measurements 
Empirical CDF of Azimuthal Angle 


Empirical CDF of Power Fraction Difference between 1*' and 2"* Beams 


— NYU Measurement 
| —- Location 3 

—— Location 2 
| —— Location 1 
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9G RADIO NETWORK PLANNING 
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5G Deployment Scenarios | 
¢ Ultra-dense networks for provision 
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° Up to hundred(s) of small power base stations 
in 1km? urban area 





¢ Multi-threading required to predict 
multiple base stations simultaneously 


* Coverage of radio network limited 
by signal and interference power 
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¢ Strong signal-to-noise-and-interference-ratio (SNIR) 
requirements for high data rates 


. 3.5 GHz frequency bands for area-wide services 
and the 26/28 GHz bands for capacity hotspots 





* Network planning allows to simulate the coverage 
before the deployment = 5G deployment strategies 
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зр DOM | ΝΑΝΤ ΡΑΤΗ MOD El Typical Channel Impulse Response 





e One path tends to carry the bulk of the power. 
Sufficient for propagation and coverage analysis. 






Г One path is u 
| dominates LN 


* Determination of path with full 3D approach 


(dB uV/m) 
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* Unlimited number of interactions (changes of 
orientation) 
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* Parameters of path determined (e.g. length, number of x ЧЩ e 
interactions, angles...) and used to compute path loss 
with semi-deterministic equations 








* Optional consideration of wave guiding possible (wave 


in | Full 3D approach 
guiding factor, based on reflection loss of walls) 


° Short prediction time 
e High accuracy 
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WINPROP SOFTWARE SUITE 


3D Dominant Path Model 


* Unique model combines high accuracy 
with short computation time for large scenarios 


F 
Altair 


Received Power 
[dBm] 


* Used by multiple operators and 
won benchmark from German Telekom 


* Wide frequency range (100 MHz up to 100 GHz) 


° For all types of cells (macro, micro, small, ...) with ү" 
Tx antennas above / below / on rooftop level "^ 


-50.00 
-55.00 





-60.00 
-65.00 
-/0.00 
-75.00 
-80.00 
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* For planning & optimization of radio networks or 
evaluation of different deployment strategies 
* Prediction of radio coverage 
* Received power 
e Path loss 
e Field strength 
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9G RADIO NETWORK PLANNING 


5G Deployment Scenarios Il 
e Beamforming on base station side 


№ 


Altair 





e Increase Rx power levels and SNIR 


Received Power 








for dedicated user [dBm] 
: -25.00 
* Heduce interference for others μη 
| -35.00 
4x4 antenna matrix provides Gi 
antenna gain of 16.7 dBi -45.00 
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otation EIRP) С> Е ΜΝ š psi: - E ' 4 Ῥ ] 2 = 60.00 
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e Optional beamforming on mobile station 


"X 


° Eight antenna elements > gain 13.3 dBi 





Altair 
* Blue areas mostly filled in > 
this density of base stations 
might be sufficient. 
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CONCLUSIONS 


¢ 5G will provide higher throughputs and many new applications 
e Efficient transmission schemes with massive MIMO usage 
e Additional spectrum allocations in higher frequency bands (e.g. 26 and 28 GHz) 


e 5G device and base station antenna design in FEKO 
e New challenges due to higher frequencies 
e FEKO combines optimization and dedicated solvers for arrays and electrically large structures 
e |deal solution for 5G antenna design 


• 5G radio channel and radio coverage analysis in WinProp 
« WinProp propagation models combine high accuracy (verified with measurements) with short computation time 
* Forall types of environments: urban, dense urban, suburban, rural, industrial, indoor, tunnel, stadium,... 
* Evaluation of 3D spatial channel profiles and channel statistics for massive MIMO 
¢ WinProp 3D ray tracing model correctly predicts the mm wave propagation 
e Ultra-dense networks require fast model for the efficient network planning 


UN 
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Thank You for Your Attention 


www.altairhyperworks.com/feko 
www.altairhyperworks.com/winprop 





